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ABSTRACT 

Molecular hydrogen transitions in the sub-damped Lyman a absorber at redshift Zaba — 
2.69, toward the background quasar SDSS J123714.60-1-064759.5, were analyzed in 
order to search for a possible variation of the proton-to-electron mass ratio p over 
a cosmological time-scale. The system is composed of three absorbing clouds where 
137 H 2 and HD absorption features were detected. The observations were taken with 
the Very Large Telescope/Ultraviolet and Visual Echelle Spectrograph with a signal- 
to-noise ratio of 32 per 2.5 km s“^pixel, covering the wavelengths from 356.6 to 409.5 
nm. A comprehensive fitting method was used to fit all the absorption features at 
once. Systematic effects of distortions to the wavelength calibrations were analyzed 
in detail from measurements of asteroid and ‘solar twin’ spectra, and were corrected 
for. The final constraint on the relative variation in ft between the absorber and the 
current laboratory value is Apt/ft = (—5.4 ± 6.3stat ± 4.0syst) x 10“®, consistent with 
no variation over a look-back time of 11.4 Gyrs. 

Key words: methods: data analysis - quasars: absorption lines - cosmology: obser¬ 
vations. 


1 INTRODUCTION 

The investigation of highly redshifted absorption systems 
in the line-of-sight of quasars, by means of high-resolution 
spectroscopic observations, has become an established and 
powerful method to constrain a possible variation of the 
laws of physics over cosmological time. The variation of di¬ 
mensionless fundamental constants, appearing as key build¬ 
ing blocks in the laws of physics, is targeted in a com¬ 
parison between astrophysical observation with labora¬ 
tory measurement. For example, the fine-structure constant 
a = j{4^Tr€o?ic) is probed via metal absorption search¬ 
ing for temporal (Webb et al. 1999) and spatial variations 
(Webb et al. 2011). 

The other dimensionless constant of nature deter¬ 
mining the structure of molecular matter, the proton-to- 
electron mass ratio fi = mp/rrif,, may be targeted by ob¬ 
serving absorption lines of a wide variety of molecules 
(Jansen et al. 2014). While radio astronomical observations 
of ammonia (Murphy et al. 2008; Kanekar 2011) and of 
methanol molecules (Bagdonaite et al. 2013b, a; Kanekar 
et al. 2015) have proven to be very sensitive, constraining 
at the level of ~ 10“^, there are are only two radio 
sources at extragalactic distances where these molecules are 
observed, and these lie at redshifts z < 1. 


Molecular hydrogen (H 2 ) is a target to investigate a pos¬ 
sible variation of pi at higher redshifts, where it is observed 
in larger numbers of absorption systems. H 2 has many spec¬ 
tral lines that are sensitive to a variation in pi (Thompson 
1975; Varshalovich & Levshakov 1993; Ubachs et al. 2007). 
Molecular hydrogen absorbing galaxies often display up to 
100 H 2 spectral lines, which helps in improving the statisti¬ 
cal basis of measurements made with this technique. Most 
importantly, H 2 absorbing systems can be observed, using 
ground-based telescopes, at redshifts z > 2, even extending 
to z = 4.22 (Bagdonaite et al. 2015), where the Lyman and 
Werner bands are redshifted into the optical band. 

So far eight H 2 absorbing systems at z > 2 have been 
analyzed for /r-variation, most notably and most accurately 
the J2123-005 system at z=2.05 (Malec et al. 2010; van 
Weerdenburg et al. 2011), the HE0027-1836 (Rahmani et al. 
2013) and the Q2348-011 system (Bagdonaite et al. 2012) 
both at 2=2.42, the Q0405-443 system at 2=2.59 (Rein¬ 
hold et al. 2006; King et al. 2008; Thompson et al. 2009), 
the Q0642-504 system at 2=2.66 (Bagdonaite et al. 2014; 
Albornoz Vasquez et al. 2014), the Q0528-250 system at 
2=2.81 (Reinhold et al. 2006; King et al. 2008, 2011), the 
Q0347-383 system at 2=3.02 (Reinhold et al. 2006; King 
et al. 2008; Thompson et al. 2009; Wendt & Molaro 2011, 
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2012), and the Q1443+272 system at 2=4.22 (Bagdonaite 
et al. 2015). A general conclusion from these studies is 
that the proton-to-electron mass ratio is constrained to 
|A/r//r| < 10“® for redshifts in the range 2=2-4.2, corre¬ 
sponding to look-back times of up to 90 per cent of the 
age of the Universe. A review of the results of the analyses 
of H 2 absorbing systems has been given by Ubachs et al. 
(2011) and by Wendt (2014). 

In the present study H 2 absorption at 2=2.69 
is investigated in the line of sight toward quasar 
SDSS J123714.60+064759.5, hereafter J1237+0647, for fur¬ 
ther constraining a variation of the proton-to-electron mass 
ratio. The importance of this system is that carbon monox¬ 
ide (CO) absorption in the optical domain is observed along¬ 
side with H 2 (Noterdaeme et al. 2010). This may be the 
basis for a simultaneous ^-variation analysis using H 2 and 
CO (Salumbides et al. 2012). As a first step the H 2 absorp¬ 
tion in this system is analyzed for constraining A/r/^, results 
of which are presented in this paper, and this will serve as 
the basis for an analysis of the CO lines in a later paper. 


2 DATA 

The analysis presented in this work is based on observa¬ 
tions made under four different observing programs, all car¬ 
ried out using the Ultraviolet and Visual Echelle Spectro¬ 
graph (UVES) mounted on the 8.2m Very Large Telescope 
(VLT) at Paranal, Chile (Dekker et al. 2000). UVES is 
a cross-dispersed echelle spectrograph with two arms that 
are functionally identical: one covers the wavelengths in the 
range 300-500 nm (Blue) and the other covers the range 
420-1100 nm (Red). Two dichroic beam splitters can be 
used to work in parallel with the two arms, centering them 
at different wavelengths. The details of the observational 
campaigns are presented in Table 1. 

The raw 2D data were reduced following the procedure 
used by Bagdonaite et al. (2014). The Common Pipeline 
Language version of the UVES pipeline was used to bias 
correct and flat field the exposures and then to extract 
the flux. The ThAr lamp exposures were used for wave¬ 
length calibration; the ThAr flux was extracted by using 
the same object profile weights as in a corresponding quasar 
exposure. After the standard reduction, the custom software 
UVES_popler^ was used to combine the extracted echelle 
orders into a single ID spectrum which was then manually 
inspected and cleaned from bad pixels and other spectral 
artifacts, and the continuum was fitted with low-order poly¬ 
nomials. 

The ‘final spectrum’, obtained by combining all 
the exposures together, covers the wavelengths from 
329.0 to 960.0 nm with a signal-to-noise ratio (S/N) of 
26 per 2.5kms“^per pixel at ~400 nm in the continuum. 


2.1 Data from 2013 and 2014 

The quasar J1237+0647 was observed in visitor mode in 
May 2013, program 091.A-0124(A), and in service mode in 
the period March-June 2014, program 093.A-0373(A). The 


total integration time amounted to 11.5 hrs. Each exposure 
was nominally 4800 s in duration and was expected to de¬ 
liver a S/N of 13 per 2.5kms“^per pixel at 400 nm and a 
resolving power R of ~40 000, with a seeing of 0.8 arcsec, 
a slit of 1.0 arcsec, an airmass of 1.2 and 2x2 binning. Due 
to bad weather conditions, four exposures were shorter than 
scheduled and, therefore, the predicted S/N level could not 
be reached in those cases. 

Each of the quasar exposures had an ‘attached’ stan¬ 
dard ThAr calibration, which is a lamp exposure taken im¬ 
mediately after the science exposure leaving all the instru¬ 
ment parameters unchanged, and was immediately followed 
by a ‘supercalibration’ exposure without allowing any grat¬ 
ing reset. The supercalibration is a method to quantify any 
long-range wavelength distortion and it involves observing 
objects with a solar-like spectrum, like asteroids or ‘solar 
twin’ stars (see Section 5). 

2.2 Data from 2009 

Part of the total dataset used in this work was retrieved 
from the ESO data archive^. This was the case for obser¬ 
vations from programs 082.A-0544(A) and 083.A-0454(A), 
which were carried out in March-April 2009 and were re¬ 
ported by Noterdaeme et al. (2010). The main difference 
with the most recent observations of J1237+0647 is that 
these exposures only had the regular ‘non-attached’ ThAr 
calibrations taken at the end of each night and dedicated ob¬ 
servations of supercalibration targets were not carried out. 
These observations of J1237+0647 from 2009 added 8.5 hrs 
of integration on the target. 


3 METHOD 
3.1 Theory 

Thompson (1975) originally proposed that a possible 
temporal variation of the proton-to-electron mass ratio 
/r = mp/rrie can be detected using absorption spectra of 
molecular hydrogen observed at high redshift. The observed 
wavelength Af of the i-th transition will show a shift given 
by: 


A/ = A°(l + 2abs)(l + Ai^), 


( 1 ) 


where A® is the rest wavelength, 2abs is the redshift where the 
absorption occurs, A/r = — /ro is the difference between 

the proton-to-electron mass ratio in the system observed and 
the value measured in the laboratory, and Ki is the sensitiv¬ 
ity coefficient which, for a given transition, determines the 
shifting power and sign due to varying /r. These coefficients 
are defined as 

A, = ^^. (2) 

d In/i 

The sensitivity coefficients for the H 2 molecule used in this 
analysis were calculated within a semi-empirical framework 
by Ubachs et al. (2007), including effects beyond the Born- 
Oppenheimer approximation. 


http: //astronomy. swin. edu. au/-iimiurphy/UVES_popler/ 


http://archive.eso.org/eso/eso_archive_main.html 
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Table 1. Details of the J1237+0647 quasar observations with UVES/VLT included in the present analysis. The total integration 
time, summed over 17 exposures, is 19.9 hrs. Most of the listed exposures were followed by attached ThAr calibrations and additional 
‘supercalibrations’ as indicated (see Section 5). The rest of the data were calibrated using the regular ThAr exposures taken at the end 
of the night. The CCD binning was 2x2 and the slit width was 1.0 arcsec for all frames. Data collected under programs 082.A-0544(A) 
and 083.A-0454(A) were retrieved from the ESO archive. 


Program ID 

Date 

Integration time [s] 

Central A [nm] 

ThAr 

Supercalibration 

082.A-0544(A) 

27-03-2009 

5400 

390 + 564 

Regular 

No 


27-03-2009 

5400 


Regular 

No 


29-03-2009 

5400 


Regular 

No 


29-03-2009 

5400 


Regular 

No 

083.A-0454(A) 

27-04-2009 

4500 

390 + 775 

Regular 

No 


27-04-2009 

4500 


Regular 

No 

091.A-0124(A) 

15-05-2013 

4800 

390 + 580 

Attached 

Yes 


15-05-2013 

1727 


Attached 

No 


15-05-2013 

4800 


Attached 

Yes 


16-05-2013 

4800 


Attached 

Yes 


16-05-2013 

2025 


Attached 

No 

093.A-0373(A) 

23-03-2014 

4800 

390 + 580 

Attached 

Yes 


03-04-2014 

4800 


Attached 

Yes 


28-05-2014 

4800 


Attached 

Yes 


29-05-2014 

826 


Attached 

No 


31-05-2014 

4800 


Attached 

Yes 


02-06-2014 

2800 


Attached 

No 


The H 2 laboratory wavelengths, A®, used as a reference 
in this analysis were measured by Salumbides et al. (2008) 
with a fractional wavelength accuracy AA/A ~ 5 x 10“®, and 
they can be considered exact in comparison with the uncer¬ 
tainties of the lines in the quasar spectrum. The complete 
list of laboratory data required for fitting the H 2 absorption 
lines observed in this work, including rest wavelengths A°, 
oscillator strengths fi, damping coefficients Ti and the sen¬ 
sitivity coefficients Ki, was compiled by Malec et al. (2010). 

3.2 Fitting Method 

A comprehensive fitting method 

(Malec et al. 2010; King et al. 2011; Bagdonaite et al. 2014), 
which involves a simultaneous treatment of all the lines, 
was used to model the H 2 spectrum. Furthermore this 
technique allows to tie part of the fitting parameters 
together, resulting in a smaller number of free parameters. 
Simultaneous fitting and parameter tying allows to include 
also those H 2 transitions that are overlapped by metal or 
H I lines from the Lyman a forest which can be modelled 
at the same time. 

The program used is the non-linear least-squares Voigt 
profile fitting program vpfit^, developed specifically for 
quasar spectra analysis. A Voigt profile is obtained from 
the convolution of a Lorentzian profile, which describes the 
natural line broadening of the transitions and is specific for 
each molecular transition considered, and a Gaussian profile, 
which reflects the physical conditions within the absorbing 
clouds by describing the velocities, both thermal and tur¬ 
bulent, of the molecules. These two profiles are eventually 
convolved with a model for the instrumental profile, which 
is assumed to be Gaussian. 

In VPFIT, the profile of each velocity component, of each 
^ http://www.ast.cam.ac.uk/-rfc/vpfit.html 


transition, is described by three free parameters: the column 
density N, the redshift 2; at which the absorption occurs, and 
the Doppler width b. These parameters are used in addition 
to the molecular physics input for every i-th transition, like 
/ii whose product with N gives the optical depth of the 
absorption line, and T^, which defines the lifetime broad¬ 
ening. Transitions probing the same rotational ground state 
share the same population, and as a consequence they are all 
described using the column density Nj. An important un¬ 
derlying assumption is that absorbing features at the same 
redshift originate from the same cloud, hence they share 
certain parameters related to the physical conditions of the 
cloud. In particular, they are assumed to share the redshift 2 ; 
and the width h, which reflects the turbulent motion and the 
kinetic temperature of the gas in the cloud. This is achieved 
by tying in vpfit the 2; and b parameters of all the ab¬ 
sorption lines associated with the gas in a single absorbing 
cloud. 

The fitting program finds the best match between the 
data and the model iteratively by minimizing the param¬ 
eter. The initial values of the free parameters are inserted by 
the user. Then, in each iteration, vpfit changes the values of 
the free parameters of the model and checks the changes in 
relative yf, reporting convergence once a stopping criterion, 
which is user-defined, is met. 

Since introducing too many free parameters may lead 
to an overfit of the spectrum, models with a different num¬ 
ber of velocity components per H 2 transition were compared 
using the yf and the Akaike Information Criterion {AICC, 
Akaike 1974) parameter. In particular, the latter parameter 
is defined as 

AICC = x^ + 2p+ (3) 

n — p — 1 

where p is the number of free parameters and n is the num¬ 
ber of spectral points included in the fit. AAICC> 5 and 
AAICC> 10 are considered to be a strong and a very strong 
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evidence respectively that the model with the lower AlCC 
is statistically preferred. 

Multiple fits, with a different number of velocity compo¬ 
nents each, were performed in order to develop a robust ab¬ 
sorption model. Only after the results of the fitting process 
are stable, the possible variation of the proton-to-electron 
mass ratio expressed by the variable A/i//i is invoked as an 
extra free parameter in the fit besides Nj, z and b. This 
fourth parameter is not introduced earlier during the fitting 
process to avoid that any over/underfitting can be compen¬ 
sated by an artificial variation of /r. The program derives 
the value of Afi/jj, according to Eq. (1), therewith also us¬ 
ing the sensitivity coefficients of Eq. (2). The consistency of 
the derived Afi/fi constraint is then tested for sensitivity to 
a number of assumptions made in the fitting process. The 
results will be discussed in detail in Section 4.2. 


4 ANALYSIS 

Quasar J1237+0647 is located at redshift 2=2.78, thereby 
defining the extent of the Lyman a forest, which comprises 
multiple series of neutral hydrogen transitions arising from 
the intergalactic medium at all redshifts below the quasar 
redshift. Toward the quasar there is one major absorp¬ 
tion system at redshift 2 = 2.69, which is a sub-damped 
Lyman a system (DLA) with a neutral hydrogen column 
density of log A = 20.0 ± 0.15 cm“^, that features many ab¬ 
sorption lines from atomic, ionized and molecular species 
(Noterdaeme et al. 2010). 

Molecular hydrogen absorption lines associated with the 
DLA are spread in the range from 352.1 to 409.5 nm. H 2 is 
found in three absorption features at redshift 2 = 2 . 688001 , 
2 =2.68868 and 2=2.68955; their velocities relative to the 
redshift of the strongest feature at 2=2.68955 are respec¬ 
tively 1 )=—125kms“^, hereafter Zl, n=—72kms“^, here¬ 
after Z2, and j;=0kms“^, hereafter Z3. The presence of 
these three features close to each other results in the line 
profiles shown in Fig. 1. The presence of clearly separated 
velocity features is not uncommon and has been observed 
in various systems in the past. The number of such compo¬ 
nents can vary from one up to the seven distinct features 
observed in the absorbing system toward quasar Q2348-011 
(Petitjean et al. 2006; Bagdonaite et al. 2012). 

Z3 has the highest column density for every J- 
level, resulting in heavily saturated absorption features 
for low J-levels with 1. Furthermore, Z3 shows also 
deuterated molecular hydrogen (HD) absorption, where 
only the J=0 level was detected (Noterdaeme et al. 2010). 
The HD molecule is sensitive to a variation of /i 
(Ivanov et al. 2008, 2010), and its detected transitions were 
included in the present ^-variation analysis. 

4.1 Creating an absorption model 

4-1.1 Selecting and fitting spectral regions 

Given the broad absorption profile, spanning ~ 150kms“^, 
overlaps between profiles of different molecular hydrogen 
transitions are common as well as overlaps with the Hi 
transitions from the Lyman a forest. Within the compre¬ 
hensive fitting method, all kinds of overlaps can be han- 



Relative velocity [km/s] 

Figure 1. Typical line profile for the three absorption fea¬ 
tures toward J1237+0647. The transition plotted is L1P(3) at 
Aq = 109.978 nm. The vertical dashed lines (red) indicate the po¬ 
sitions of the three H 2 spectral features fitted in this study. The 
velocity frame is centered on the strongest feature Z3. 

died. However, no relevant information was gained in the 
case of a complete overlap with a saturated Hi line, hence 
such absorption features were not included in this study. 
Less frequently, some H 2 features partially overlap with nar¬ 
row metal transitions. However, these lines were not ex¬ 
cluded from the dataset considered in this work. A large 
number of metal absorption features were identified by 
Noterdaeme et al. (2010). 

All the 137 potentially useful H 2 and HD transi¬ 
tions were contained in 66 spectral regions, in the range 
351.0-413.5 nm, which were selected trying to avoid, when¬ 
ever possible. Hi transitions and too-heavily saturated low 
J-level H 2 absorption features in Z3. The useful H 2 transi¬ 
tions, probing six different rotational quantum states with 
J=0-5, are listed in Table 2. Absorption features originating 
at the same redshift had both the 2 and the b parameters 
tied together, while the column densities N j were allowed to 
vary independently from each other. Z3 shows saturated H 2 
lines, particularly in the low J-levels, but there is an appar¬ 
ent evidence of a partial coverage with a non-zero residual 
flux at the base of its absorption features. To take into ac¬ 
count this effect, a zero-level correction was included as a 
free parameter in the fit for each H 2 region. 

In order to build a robust model, the neutral hydrogen 
absorption lines occurring in the Lyman a forest and in¬ 
cluded in the selected regions must be accounted for. The 
Hi transitions were modelled by assigning to each of them 
a set of free parameters {N, 2, b) in vpfit. Initial guesses of 
the free parameters were user-provided for each H 2 , HD and 
Hi transition. 

4-1.2 Extra velocity components 

The H 2 absorption spectrum toward J1237+0647 displays 
three clearly distinct features in velocity space. In order 
to fit the spectrum as accurately as possible, it is impera¬ 
tive to consider possible underlying substructures in velocity 
space. Physically this may relate to the existence of closely 
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Table 2. List of the molecular transitions used in this analysis. 


J-level 

Lyman transitions 

Werner transitions 

^transitions 

Ha J = 0 

LOR(O), L2R(0), L3R(0), L4R(0), L6R(0), L7R(0), 
L8R(0), L10R{0) 

W0R{0), WIR(O) 

10 

H 2 J = 1 

LOP(l), LOR(l), L1P{1), LlR(l), L2P(1), L2R(1), L3P(1), 
L7P(1), L9R(1), LlOR(l), L12R(1) 

W0R{1), WlR(l) 

13 

Ha J = 2 

L0P(2), L0R(2), L1P(2), L1R(2), L2P(2), L2R(2), L3P(2), 
L4P(2), L4R(2), L5P(2), L6P(2), L7P(2), L7R(2), L8P(2), 
L9P(2), L10P(2), L10R(2), L11P(2), L11R(2), L13P(2), 
L13R(2) 

W0P(2), W0Q(2), W1Q(2), W2P(2), 
W2Q(2), W2R(2) 

27 

Ha J = 3 

L0P(3), L0R(3), L1P(3), L1R(3), L2R(3), L3P(3), L3R(3), 
L5P(3), L5R(3), L6P(3), L6R(3), L7P(3), L7R(3), L8P(3), 
L8R(3), L9R(3), L10P(3), L12R(3), L13P(3i L13R(3), 
L14P(3) 

W0P(3), W0Q(3), W1R(3), W2P(3), 
W2Q(3), W2R(3), W3P(3) 

28 

Ha J = 4 

L0P(4), L0R(4), L1P(4), L1R(4), L2R(4), L3P(4), L3R(4), 
L4P(4), L4R(4), L5R(4), L6P(4), L7P(4), L7R(4), L8P(4), 
L9P(4), L9R(4), L11P(4), L11R(4), L12R(4), L13P(4), 
L14P(4) 

W0P(4), W0Q(4), W0R(4), W1P(4), 
W2Q(4), W2R(4), W3Q(4) 

28 

Ha J = 5 

L0R(5), L1P(5), L2P(5), L2R(5), L3P(5), L3R(5), L4P(5), 
L5P(5), L5R(5), L7P(5), L8P(5), L9P(5), L10P(5), 
L11R(5), L12R(5), L13P(5), L14P(5), L15P(5) 

W0P{5), W0R(5), W1Q(5), W2Q(5), 
W2R{5), W3P(5) 

24 

HD J = 0 

L3R0, L5R0, L8R0, L12R0, L15R0 

WIRO, W3R0 

7 

Total number of lines 


137 


separated distinct clouds, or inhomogeneities in the clouds 
in the line-of-sight. To investigate such underlying velocity 
structure, the fitting residuals from the 66 selected regions 
were normalized, shifted to a common velocity scale and 
averaged together, creating a composite residual spectrum 
(CRS, Malec et al. 2010). Only the H 2 absorption features 
that were not overlapped by Hi features were included in 
the CRS, in order to avoid that the residuals of the Hi 
model could introduce any spurious effect. The CRS, by 
combining the residuals of many transitions, highlighted any 
over/underfitted structure in the absorption features. 


Z1 could be described using only one velocity compo¬ 
nent (VC), while indications of underfitting were observed 
in the CRS for absorption features Z2 and Z3, hence they 
were modelled by adding extra VCs. They were added to the 
model until vpfit started to reject them. Then the models 
with a different number of VCs were compared using the 
and the AICC parameters, in order to decide which one 
better fitted the data. The CRS for the three absorption 
features Zl, Z2 and Z3 are shown in Figs. 2-4. Z2 could 
be modelled using two VCs, labelled as ‘a’ and ‘b’ in Fig. 
3, and Z3 was modelled using three VCs, labelled as ‘a’, 
‘b’ and ‘c’ in Fig. 4. Parameters for the model are summa¬ 
rized in Table 3 and the VCs are shown in Fig. 5. Previ¬ 
ous studies found no strong evidences for a systematic shift 
in Afi/^ due to a segregation in different rotational levels 
(Malec et al. 2010; King et al. 2011), hence it was assumed 
that the velocity structure is the same among all the J-levels 
in each absorption feature. The complete absorption model 
is shown in Section A. 



Figure 2. Top panel: normalized composite residual spectrum for 
Zl formed from 30 H 2 absorption features. The dashed lines rep¬ 
resent the ±1(T boundaries. Bottom panel: an example transition 
shown on the same velocity scale. The velocity scale is centred at 
z = 2.688009. The vertical dashed line shows the position of the 
velocity component. 

4.2 Constraining A^//r 

Once an optimal model was found, an extra free param¬ 
eter, corresponding to A/r//i, was added to the H 2 tran¬ 
sitions in VPFIT. All the transitions were included in the 
ht, without any distinction between their bands, their J- 
levels or absorbing clouds where they arose. The model re¬ 
turns = (—1-1 ± 6.3stat) 10“®. The statistical error 

reported here was derived from the appropriate diagonal 
term of the final covariance matrix for the fit. Given the ab¬ 
sorption model, it represents only the uncertainty in A/r//r 
stemming from the photon statistics of the quasar spectrum, 
i.e. its signal-to-noise ratio. However, to evaluate the robust- 
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Table 3. Column densities Nj^ redshift z and Doppler parameters h of the H 2 transitions in the three clearly distinguishable absorption 
features and in their underlying velocity components a, b, c. Levels with J = 4 and 5 are too weak to be detected in Z1 and Z2. Levels 
with J = 0 and 1 are heavily saturated in Z3 and they have been discarded from the model. 



Zl 

Z2 



Z3 


J-levels 

logAj [cm-2] 

logAj [cm 2] 


logAj [cm-2] 




a 

b 

a 

b 

c 

J = 0 

15.74 ± 0.11 

15.44 ± 0.11 

13.78 ±0.17 


Discarded 


J = 1 

16.29 ±0.08 

16.23 ±0.11 

14.56 ±0.10 


Discarded 


J = 2 

15.35 ±0.04 

15.27 ±0.07 

13.90 ±0.15 

15.67 ± 0.10 

17.64 ± 0.10 

15.26 ±0.16 

J = 3 

15.21 ± 0.02 

14.88 ±0.05 

14.31 ± 0.08 

15.66 ±0.11 

18.02 ± 0.06 

14.64 ± 0.28 

J = 4 


N/D 


14.79 ±0.07 

15.34 ±0.04 

13.93 ±0.07 

J = 5 


N/D 


14.50 ±0.07 

14.79 ± 0.04 

13.33 ±0.17 

HD J = 0 


N/D 


N/D 

14.36 ± 0.03 

N/D 

z 

2.68801(06) 

2.68866(09) 

2.68872(82) 

2.68950(96) 

2.68955(14) 

2.68972(52) 

b [kms“^] 

2.94 ± 0.08 

2.32 ±0.13 

6.99 ± 0.65 

11.36 ±0.48 

5.78 ±0.11 

2.54 ± 0.48 




Relative velocity [km/s] Relative velocity [km/s] 


Figure 3. Top panel: normalized composite residual spectrum 
for Z2 formed from 33 H 2 absorption features obtained from the 
2 VCs model. The dashed lines represent the dilcr boundaries. 
Bottom panel: an example transition shown on the same veloc¬ 
ity scale. The velocity scale is centred ot z = 2.688664, on the 
strongest velocity component. The vertical dashed lines show the 
positions of the two velocity components at u = Okms“^(a) and 
-|-7kms“^(b). 


Figure 4. Top panel: normalized composite residual spectrum 
for Z3 formed from 33 H 2 absorption features obtained from the 
3 VCs model. The dashed lines represent the dilcr boundaries. 
Bottom panel: an example transition shown on the same veloc¬ 
ity scale. The velocity scale is centred at 2 : = 2.689555, on the 
strongest velocity component. The vertical dashed lines show the 
positions of the three velocity components at u = —4kms“^(a), 
Okms“^(b) and -|-13 km(c). 


ness of the result, hereafter the statistical value of a 

variety of consistency tests and sources of systematic error 
which contribute to the overall uncertainty budget were ex¬ 
plored, results of which are shown in Fig. 6 and discussed 
below. 


4 . 2.1 Isolating single absorbing features 

Each of the absorbing features Zl, Z2 and Z3 can 
be investigated separately. The constraints on A/i/^ 
returned are = (11.4±10.2stat) X IO-® for 

Zl, A/r//r|z 2 = (-17.7±14.0stat)xl0-® for Z2 and 
A/r//i|z 3 = (—6.1 ± 14.8stat) X 10 ® for Z3. The resulting 
constraints are in agreement with each other and with the 
statistical value for the entire dataset within Ict boundaries. 


4-2.2 Isolating Lyman and Wemer transitions 

All the H 2 transitions considered in this analysis belong ei¬ 
ther to the Lyman- or to the Werner band systems. Any 
long-range distortion of the wavelength scale will be nearly 
degenerate with a shift of the molecular hydrogen lines due 
to a non-zero value of Afi/fi. This degeneracy is partly bro¬ 
ken by fitting together H 2 transitions from both Lyman and 
Werner bands. The effect of such a distortion can be inves¬ 
tigated by fitting separately transitions from the two band 
systems. 

The values returned in this way are 
A/t/^Il = (—0.6 ± 6.4stat) X 10“® for Lyman transi¬ 
tions and Afj,/g-\yj- = (—10.3 ± 33.9stat) x 10“® for Werner 
transitions. The larger error of the value returned by 
Werner transitions is due to three effects. Firstly, there are 
fewer Werner transitions than Lyman transitions in the 
spectrum, with a ratio between the former and the latter 
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Figure 5. Example transitions for every J-level. Each velocity component is plotted against the fitted continuum. Hi contributions 
to the line shapes are not highlighted. Top panel: J-levels 0-4 for Z1 (blue), the velocity scale is centred at ^=2.688009. Middle panel: 
J-levels 0-4 for Z2. The velocity component at z=2.688664, Z2a, is depicted with the (blue) solid line and the one at z=2.688717, Z2b, 
with the (red) dashed line. The velocity scale is centred on Z2a. Bottom panel: J-levels 2-5 for Z3. The three velocity components at 
redshifts z=2.689498, Z3a, 2.689555, Z3b, and 2.689723, Z3c, are shown with the dashed (red), solid (blue) and dotted (green) line 
respectively. The velocity scale is centred on Z3b. 


of ~1:4. Then the spread of the sensitivity coefficients is 
larger for Lyman transitions, AK]^ = 0.092, than the one 
of Werner transitions, AKt^ = 0.054. Finally, the Werner 
transitions fall in the bluest part of the spectrum, where 
the S/N is much lower than in the region where Lyman 
transitions are. Such differences lead the statistical value of 
An/fi to be dominated by the Lyman transitions. 


4-2.3 Isolating cold and warm states 

Using VPFIT it is possible to calculate a different value of 
A/a/fi for transitions originating in each J-level, in order to 
obtain the relative contributions of different levels to the 
statistical value. This is particularly useful to investigate 
the impact of possible temperature inhomogeneities in the 
absorbing clouds. Due to the para-ortho distribution, the 
J= 1 rotational state is significantly populated even at the 
lowest temperatures (Ubachs et al. 2007). As a consequence, 
the H 2 transitions can be divided into two sets: rotational 
states with J=0,1 (cold) and with J)S 2 (warm). 

The constraints returned by the model are 
= (—15.1 ± 11.5stat) 10“® for the cold states 

and A/r//r|warm = (4.4 ± 7.1stat) X for the warm 

states. The two values match within their uncertainties, 
which means that any temperature inhomogeneity in the 


absorbing clouds does not have a significant impact on the 
p-variation. 

Furthermore, the slightly larger error derived from the 
cold states reflects the fact that levels with J^ 3 are satu¬ 
rated in Z3 and, in particular, cold states are heavily satu¬ 
rated. As a consequence of the saturation, there are fewer 
constraining pixels contributing to the fitting process and 
there are fewer transitions in the cold states dataset. 


4-2-4 Separating exposures from 2009 

The dataset used in this work can be divided between the 
exposures taken in 2013 and 2014, which had an attached 
ThAr calibration and a supercalibration, and the exposures 
taken in 2009, retrieved from the ESO archive, which only 
had regular non-attached ThAr calibrations. In order to see 
how the absence of the individual calibrations affects the 
statistical value, two values for Ap,/p, were derived from the 
two subsets. 

The constraints obtained in this way are 
A/r /^|2009 = (“9-1 i 7.6stat) X 10“® for the exposures 
taken in 2009 and A^//r|2014 = (14.3 ± 10.3stat) x 10~® for 
the exposures taken in 2013 and 2014. The two values do 
not agree within their combined la statistical uncertainties, 
although they do within 1.5a. However, this comparison 
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Figure 6. constraints obtained from various tests (trian¬ 

gles) compared with the statistical value from the entire dataset 
(circle). The solid (green) horizontal line indicates the value of 
the statistical constraint and the two dashed lines represent its 
1(7 statistical uncertainties. The square represents the A/j/^ value 
obtained only from the exposures taken in 2014 after applying to 
each of them a long-range distortion correction, as discussed in 
Section 5. The diamond (red) represents the fiducial value ob¬ 
tained for the entire dataset (2009-1-2013-1-2014) after applying 
the distortion corrections. 

does not take into account the long-range distortion efTect, 
which will be discussed in the next section. 


5 SUPERCALIBRATIONS 

In order to constrain a variation of fi, it is crucial to 
accurately wavelength calibrate the quasar observations. 
Attempts to perform advanced tests of the effective ac¬ 
curacy of UVES were performed by Molaro et al. (2008) 
by comparing a UVES-recorded asteroid spectrum with 
a highly accurate solar reference spectrum. These mea¬ 
surements revealed zero offsets up to ~ 50 ms“^, likely 
due to a non-uniform slit illumination, but found no 
evidence for a wavelength-dependent velocity shift in 
the spectrum. Similar studies, that are now referred to 
as ‘supercalibrations’, showed that spectra taken with 
UVES and HIRES spectrographs suffer from intra-order 
calibration shifts on scales of single echelle orders but 
they did not uncover any evidence for a long-range 
distortion in the spectra on the scale of the spectro¬ 
graphs arms (Whitmore et al. 2010; Griest et al. 2010). 
Rahmani et al. (2013) detected for the first time 
a long-range distortion in some UVES spectra up 
to ~ 400 ms“^ on a scale of ~ 600 A. More recently, 
Whitmore & Murphy (2015) found that significant wave¬ 
length calibration distortions, both intra-order and long- 
range, are ubiquitous across UVES and HIRES history, 
most likely due to a different light path of the astronomical 
object and the ThAr lamp within the instrument. 

The presence of long-range distortions between the at¬ 
tached ThAr calibration and the main target spectrum ve¬ 
locity scales introduces a wavelength-dependent shift that 
mimics a non-zero A/i//i. This happens because the sen¬ 


sitivity coefficients Ki decrease with increasing wavelength 
(see Fig. 1 of Bagdonaite et al. (2014)) and a long-range 
distortion of the wavelength scale would produce a strong 
systematic effect which will be nearly degenerate with the 
effect produced by a variation of /r. In principle, this degen¬ 
eration can be broken by fitting Lyman and Werner tran¬ 
sitions together; however, the paucity of Werner transitions 
and the low S/N in the spectral region where they fall pre¬ 
vent the degeneracy to be broken in this way, as illustrated 
by Malec et al. (2010) and discussed in Section 4.2.2. Over 
most of UVES’s history, the sign of measured long-range 
distortions was such that they spuriously pushed to 

more positive values (Whitmore & Murphy 2015). 


5.1 Method 

The first supercalibration performed with UVES con¬ 
sisted of a comparison between asteroid spectra and 
a Fourier Transform Spectrometer (FTS) solar atlas 
(Molaro et al. 2008). Rahmani et al. (2013) cross-correlated 
asteroid spectra with a laboratory solar spectrum and, more 
recently, Whitmore Murphy (2015) improved the tech¬ 
nique by forward modelling the FTS spectra to match ob¬ 
served asteroid spectra, thereby allowing information to be 
derived from much shorter spectral ranges, and by showing 
that the same technique can be applied on solar twin stars 
as well. Their method is used in this study to supercalibrate 
the quasar spectrum. 

The supercalibration process consists in the comparison 
of a ThAr calibrated spectrum with a reference spectrum 
from a FTS, whose frequency scale is expected to be much 
more accurate than the UVES one, hence it is considered 
‘absolutely calibrated’ for the purpose of this study. The ref¬ 
erence used is the FTS vacuum solar spectrum from Chance 
& Kurucz (2010)"^. The targets for the supercalibrations 
were asteroids, which reflect the solar light and hence have 
the same spectrum as the Sun, and ‘solar twin’ stars, i.e. 
objects whose spectra are almost identical to the solar one 
(Melendez et al. 2009; Datson et al. 2014). Since the super¬ 
calibration targets are astronomical objects, no changes are 
required in the focus or in the slit alignment. 

The two spectra were divided into small regions 
of ~ 500kms“^ each, which corresponds to ~ 8 A. Their 
small size allows for sampling each echelle order with « 10 
regions, and for identifying any relative velocity distor¬ 
tion down to ~ 30 ms“^ (Whitmore & Murphy 2015). The 
spectral regions of the two spectra were compared using 
a x^-minimization technique looking for any wavelength- 
dependent calibration distortions across each echelle order 
and across the wavelength range covered by all the echelle 
orders collectively. The overall velocity distortion was ob¬ 
tained by a linear fit of all the spectral regions and it was 
used to distortion-correct the corresponding quasar expo¬ 
sure. 


Available at http://kurucz.harvard.edu/sun/ 

irradiance2005/irradthu.dat 
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5.2 Supercalibration data 

5.2.1 Data from 2013 and 2014 

The exposures taken in 2013 and in 2014 were recorded with 
an attached ThAr wavelength calibration. Since the prob¬ 
lems of long-range distortions and their effect on /r-variation 
analysis from H 2 spectra had been well identihed and 
reported (Rahmani et al. 2013; Bagdonaite et al. 2014), su¬ 
percalibration exposures targeting asteroids and solar twins 
were purposely recorded immediately after the science expo¬ 
sures with their attached ThAr calibrations. In view of their 
accessibility during observations of J1237-I-0647, supercali¬ 
bration spectra of the Eunomia asteroid, as well a solar twin 
stars HD 76440, HD 147513, HD 097356 and HD 117860, 
were used. Only some exposures, which had to be aborted 
due to bad weather conditions, were not supported by a su¬ 
percalibration. 

For each supercalibration spectrum, its long-range 
wavelength distortion was determined in terms of a slope 
parameter, which was found to range between 300 and 
725 ms”'^ per 1000 A. In particular, during each observ¬ 
ing run the values of the distortion slopes showed an ap¬ 
parent stability, with a spread in the slope values of, at 
most, ~200 ms“^ per 1000 A in May 2014. These val¬ 
ues, and in particular the positive sign, are commensu¬ 
rate with values found previously in such measurements 
(Rahmani et al. 2013; Whitmore & Murphy 2015). Results 
of the present analysis are plotted in Fig. 7 and numer¬ 
ical values listed in Table 4. Thereupon each individual 
J1237-I-0647 science exposure was corrected for its long- 
range distortion by introducing a counter-distortion based 
on these slopes, after which the exposures were rebinned 
and combined using UVES_popler. Exposures without a 
directly attached supercalibration were corrected using the 
average value of all the distortion slopes from the supercali¬ 
brations taken during the same night, in view of the stability 
discussed above. 

5.2.2 Data from 2009 

The data for J1237-I-0647 taken in 2009, retrieved from 
the ESO archive, were not supported by an attached ThAr 
calibration, and they were therefore wavelength calibrated 
against a regular ThAr calibration recorded at the end of 
the night. No supercalibration spectra were recorded with 
the J1237+0647 science exposures. In the ESO archive some 
suitable asteroid observations were found as three expo¬ 
sures of the Ceres asteroid, recorded in March and April 
2009 within one week of the quasar exposures, under pro¬ 
gram 080.C-0881(B) (PI Dumas). These exposures did not 
share the same telescope settings as the J1237-I-0647 ob¬ 
servations. In particular they had different grating settings, 
covering only half of the H 2 window, as shown in Fig. 7, 
and they did not have an attached ThAr calibration. Un¬ 
like the quasar observations, these exposures were not taken 
with the slit perpendicular to the horizon and hence an at¬ 
mospheric dispersion corrector was used not to lose the blue 
part of the flrrx due to the atmospheric dispersion. Such 
differences can affect the light path of Ceres within the tele¬ 
scope, resulting in a different one from the quasar light path. 
Since the causes of the long-range distortions are unclear 
(Whitmore & Murphy 2015), it is not possible to estimate 


the impact of those differences on the supercalibration pro¬ 
cess. 

Analysis of the Ceres exposures yields two distortion 
slopes of of ~ 150 ms~^ per 1000 A, while a third one shows 
a negative distortion of ~ —500 ms“^ per 1000 A. The pres¬ 
ence of two values showing distortions with opposing signs 
for the same night is remarkable, in particular since these 
two Ceres spectra were calibrated using the same non- 
attached ThAr spectrum, and it prevents one to calculate a 
reliable average value of the distortion correction for these 
exposures. For these reasons, it was decided to not perform 
a counter-distortion analysis of the wavelength scale. 

Nevertheless, a magnitude of the wavelength distortion 
can be estimated using the limiting values of the distortion 
slopes, given in Table 4. Both limiting values of opposing 
signs were considered as the largest distortions that could 
have affected the UVES exposures of J1237-I-0647 in 2009. 
Therefore the full sub-spectrum formed from 2009 exposures 
was first counter-distorted using the limit values of the dis¬ 
tortion slope and then was combined with the distortion- 
corrected spectrum formed from exposures taken in 2013 
and 2014 to derive the constraints on The two val¬ 

ues derived were Afi/n\i^^ = (—4.6 ± O.lgtat) x 10“® and 
A^/^lup = (0.1 ± 6.1gtat) X 10~® for the lower and for the 
upper limit respectively. The spread in resulting values for 
A/j./jj. was interpreted in terms of a systematic uncertainty 
amounting to 2.4 x 10“®. 

5.2.3 Other sources of systematics 

The supercalibrations allow to address the issues of both the 
intra-order and the long-range distortions that are known to 
affect UVES(Molaro et al. 2008; Griest et al. 2010; Whit¬ 
more et al. 2010; Rahmani et al. 2013; Whitmore & Murphy 
2015). All the supercalibration targets used in this analysis 
are presented in Fig. 7. In each panel, the correction that 
needs to be applied to the ThAr wavelength is shown; note 
that each echelle order is characterized by several (typically 
4-9) measurements separated by ~8 A. The vertical spread 
in each echelle order reflects the magnitude of the intra-order 
distortions, while the slopes, obtained by least-squares fits to 
the mean velocity shift correction from each echelle order as 
a function of wavelength (as in Whitmore & Murphy 2015), 
represents the long-range distortions. 

The impact of the intra-order distortions is, typ¬ 
ically, not expected to be dominant in the systematic 
error budget on the fiducial value of Aja/fi, because the 
H 2 and HD transitions are spread across multiple orders 
(Malec et al. 2010; King et al. 2011; Bagdonaite et al. 2014) 
Nevertheless, the supercalibrations taken in 2014 show 
intra-order distortions that are ~ three times larger than 
the distortions in the exposures taken in 2009 and 2013. 
An estimation of the impact of the intra-order distor¬ 
tions on the systematic error on Afj,/yL is made using 
5(A/i//i) = [{Av/c)/\fN]/AKi, where Av/c is the mean 
amplitude of the intra-order distortions, N = 137 is the 
number of H 2 and HD transitions detected and AKi = 0.065 
is the spread in the sensitivity coefficients. 

The supercalibration exposures return a mean velocity 
shift of ~ 25 ms“^ for 2009 and 2013 and ~ 60 ms“^ for 
2014. The shifts translate into systematic uncertainties on 
Afi/fi of 1.2 X 10“® and 3.2 x 10“® respectively. The aver- 
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age value of 5{Afi/fj,) — 2.2 x 10“® is taken as the systematic 
uncertainty on the /r-variation due to the intra-order distor¬ 
tions. However, this estimation does not take into account 
various factors like the different sensitivity of the transi¬ 
tions, where they fall with respect to the order centre and 
the variation of the S/N across the H 2 window. As a con¬ 
sequence, this value should be considered as an upper limit 
(Malec et al. 2010). 

Long-range and intra-order distortions are not the 
only contributors to the systematic error budget. In 
addition, other effects, like the spectral redispersion, 
the presence of exposures without any ThAr attached 
calibration and the calibration residuals, contribute to 
the total systematic error budget. Previous studies ad¬ 
dressed these possible error sources, finding that the ab¬ 
sence of attached ThAr calibrations introduces an er¬ 
ror of ~ 0.7 X 10“® (Bagdonaite et al. 2014), the calibra¬ 
tion residuals have an effect of, at most, ~ 2.0 x 10“® 
(Murphy et al. 2007; Bagdonaite et al. 2014), and the spec¬ 
tral redispersion can contribute up to ~ 1.4 x 10“® for 
a pixel size of ~ 2.5kms“^ (King et al. 2011). Adding in 
quadrature all these contributions to the systematic error 
budget yields a systematic error on the fiducial value of 
of ~ 4.0 X 10"®. 

5.3 Constraining Afj./fi from the 
distortion-corrected spectrum 

After having distortion-corrected the exposures 
taken in 2013 and 2014, an updated value of 
A/r//r| 2 oi 4 = (1.2 ± ll.Sgtat) ^ 10“® was derived from 
the sub-spectrum formed from those exposures only. This 
value, which agrees within Icr with the constraint derived 
only from exposures from 2009, is represented as a square 
in Fig. 6. 

The ‘final spectrum’, obtained combining the dis¬ 
tortion corrected exposures together with the un¬ 
corrected exposures from 2009, delivers a fiducial 
constraint on the proton-to-electron mass ratio of 
A/r//r = (— 5.4 ± 6.3stat i 4-Osyst) X 10“®, which is the 
final result of the present study. 


6 CONCLUSION 

In this work, the analysis of molecular hydrogen absorption 
in the system at redshift z = 2.69, in the line of sight toward 
quasar J1237-I-0647, is presented in order to constrain a pos¬ 
sible variation of the proton-to-electron mass ratio /r. 137 H 2 
and HD transitions, found in three distinct velocity features 
associated with the DLA, some even exhibiting further ve¬ 
locity fine-structure, were analyzed. The large number of ab¬ 
sorption features considered in the dataset includes partially 
overlapped features, due to the complex velocity structure 
of the system, as well as strongly saturated lines. This was 
possible because of the comprehensive fitting method used in 
this analysis. Intra-order and long-range distortions, which 
are known to affect the LIVES spectra, were taken into ac¬ 
count by applying the supercalibration technique presented 
by Whitmore & Murphy (2015) to the exposures taken in 
2013 and 2014. For exposures taken in 2009, and retrieved 



Wavelength [A] 


Figure 7. Distortion maps of the UVES wavelength scale de¬ 
rived from asteroid and solar twin targets used for the super¬ 
calibration of science exposures of J1237-h0647. For each super¬ 
calibration exposure, the plot shows the velocity shift required 
to align the ThAr-calibrated wavelength scale with that derived 
from the FTS solar spectrum. One observes both intra-order and 
long-range distortions of the wavelength scale. Only the slopes 
of the latter, whose values are listed in Table 4 with a one-to- 
one correspondence, are the physically relevant parameters used 
in the distortion analysis. They were determined by least-squares 
fits to the mean velocity shift correction from each echelle order 
as a function of wavelength. The two dashed vertical lines show 
the window in which the H 2 absorption lines toward J1237-I-0647 
fall. The radial velocities are unimportant for the analysis of the 
slopes and have been removed here; the distortion maps have all 
been shifted to zero velocity shift at 3800 A. Top panel: three 
supercalibrations extracted from Ceres exposures taken in 2009 
and retrieved from the ESO archive. Lower panel: seven super¬ 
calibrations extracted from asteroid and solar twin observations 
performed in 2013 and 2014. Eunomia, HD76440 and HD147513 
were observed in visitor mode in May 2013 while the solar twins 
HD097356 and HD117860 were observed in service mode in 2014. 
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Table 4. Observations with VLT/UVES of the supercalibration targets used in this work. Data from program 080.C-0881{B) were 
retrieved from the ESO archive and have only the regular ThAr calibration taken at the end of the night. The slit width was 1.0 arcsec 
for all exposures. The uncertainty on the slopes is ~ 30 ms“^ per 1000 A. 


Target 

Program ID 

Date 

Grating [nm] 

Distortion slope 
[ms“^ per 1000 A] 

Ceres 

080.C-0881(B) 

23-03-2009 

346 

120 



01-04-2009 

346 

200 



01-04-2009 

346 

-510 

Eunomia 

091.A-0124(A) 

15-05-2013 

390-1-580 

300 

HD76440 


15-05-2013 

390+580 

400 

HD147513 


16-05-2013 

390+580 

370 

HD097356 

093.A-0373(A) 

23-03-2014 

390+580 

650 



03-04-2014 

390+580 

720 

HD 117860 


28-05-2014 

390+580 

460 



31-05-2014 

390+580 

260 


from the ESO archive, the impact of such distortions was es¬ 
timated using asteroid observations taken within one week 
of the quasar ones, and it was translated into a contribu¬ 
tion to the systematic error budget. The resulting value of 
a constraint on a varying proton-to-electron mass ratio is 
= (-5.4 ± 6.3stat ± 4.0spt) x IQ-®. 

This constraint can be improved using the upcom¬ 
ing new generation of high resolution spectrographs, like 
the Echelle SPectrograph for Rocky Exoplanet and Stable 
Spectroscopic Observations (ESPRESSO, Pepe et al. 2010). 
The spectral resolution achieved in this work allows to 
resolve components with Doppler broadening parameters 
b ~ 2.5kms“^, while the high resolution of the new instru¬ 
ments, roughly three times higher than used here, will allow 
to resolve components with Doppler broadening parameters 
of the order of ~ l.Okms”^. This may help further resolve 
finer structure within individual spectral features and hence 
result in a more precise A/r//r constraint from J1237-I-0647. 
Moreover, the new instrument with fiber-fed and frequency- 
comb support should not suffer from the long-range wave¬ 
length distortions that are affecting UVES, and should lead 
to smaller systematic uncertainties in /i-variation analyses. 

In Fig. 8 the result of the present work is compared 
with other constraints, taken from literature, on a varying 
H obtained from previous investigations on other H 2 ab¬ 
sorbers at medium-to-high redshifts. The weighted mean, 
obtained considering, where possible, both the statistical 
and the systematic errors, resulting from this larger set is 
Afi/jj, = (2.9 ± 1.7) X 10“®. This is consistent with no vari¬ 
ation of fj, over a cosmological time-scale at a level of 
for a look-back time of ~10.5-12.5 Gyrs. 
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Figure 8. Overview of results from investigations on a vary¬ 
ing proton-to-electron mass ratio using molecular hydrogen ab¬ 
sorbing systems. The result on J1237-I-0647 presented in this 
work is indicated with a (red) square. The circles show results 
from previous analyses on systems: a) J2123-0050 (Malec et al. 
2010; van Weerdenburg et al. 2011), b) HE0027-1836 (Rah¬ 
man! et al. 2013), c) Q2348-011 (Bagdonaite et al. 2012), d) 
Q0405—443 (King et al. 2008), e) B0642-5038 (Bagdonaite et al. 
2014), f) Q0528-250 (King et al. 2008, 2011), g) Q0347-383 (King 
et al. 2008; Wendt & Molaro 2012) and h) Q1443-I-272 (Bagdon¬ 
aite et al. 2015). Note that multiple values for the constraint on 
Afi/ji were derived from systems J2123-0050 (labelled as ‘a’), 
Q0528-250 (labelled as ‘f’) and Q0347-383 (labelled as ‘g’). To 
avoid overlaps between observations of the same system, their 
points are presented with an offset of -1-0.05 on the 2 :-scale. The 
dotted line represents the zero level, while the dashed line shows 
the weighted mean of all the Afi/fi values and the shaded area 
shows its Itr boundaries. 
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APPENDIX A: ABSORPTION MODEL 

Figures A1-A7 show the H 2 window in the spectrum of the 
absorbing system at z ^ 2.69 toward quasar J1237+0647 
along with the absorption model. 

This paper has been typeset from a T^jX/DT^?^ file 

prepared by the author. 
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Figure Al. Spectrum of quasar J1237+0647 (part 1 of 7). Regions containing the H 2 transitions used in this analysis are indicated by 
a green line of the fitted model. H 2 transitions are marked with a vertical solid line, HD transitions are marked with a dashed line and 

metal transitions are marked with a dotted line. 
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Figure A2. Spectrum of quasar J1237+0647 (part 2 of 7), continued. 
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Figure A3. Spectrum of quasar J1237+0647 (part 3 of 7), continued. 
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Figure A4. Spectrum of quasar J1237+0647 (part 4 of 7), continued. 
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Figure A5. Spectrum of quasar J1237+0647 (part 5 of 7), continued. 
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Figure A7. Spectrum of quasar J1237+0647 (part 7 of 7), continued. 
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